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ABSTRACT 

ON THE NATURAL. TRANSVERSE MODES OF A SUBMARINE 

by 

EDWARD JOHN OTTH 
and 
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SUBMITTED TO THE DEPARTMENT OF NAVAL. ARCHITECTURE 
AND MARINE ENGINEERING ON MAY 23, 1955, IN PARTIAL 
FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF 
NAVAL ENGINEER. 

A graphical method of determining normalized curvature 
patterns, deflection patterns, and frequencies for the natural 
modes of interest of a ship treated as a non-uniform, free-free, 
transversely vibrating beam is presented. Corrections for shear 
deflection and rotatory inertia are integrated into the method, 
which is based on procedures developed by Stodola and 
Schadlofsky. Integrations are performed mechanically by an 
integraph. The method is demonstrated by determining the first 
three horizontal modes of a three -eighths scale model of a 
submerged submarine. 

The method is found to be convergent, yielding acceptable 
results, and taking account of the secondary effects considered. 
Three cycles of iteration for each of the first three non-trivial 
natural modes of the model yielded frequencies of 417, 838, and 
1391 cycles per minute, respectively. Patterns obtained are tabu- 
lated in detail. 

The process developed is recommended for consideration by 
small research facilities with limited budgets. 
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Title: Associate Professor of Mechanical Engineering 
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NOMENCLATUBE 



A * total cross-section area of hull structure considered effective 
in bending and shear, (&q, inchs) 

A * axis translation of assumed or calculated deflection for ortho- 
gonality with rigid body translation 

B * rotation coefficient for axis of assumed or calculated deflection 
for orthogonality with rigid body rotation 



S’ Si* ** 



Cj * 



G — 




ac 



adjustment coefficients for correcting assumed or calculated 
deflections for orthogonality with previous modes 

constant of integration applied to shear (equal to aero) 

constant of integration applied to bending (equal to aero) 

constant of integration applied to slope 

constant of integration applied to deflection 

6 

Young’s modulus, 30 x 10 psi 

6 

shear modulus, 11. 5 ac 10 psi 

moment of inertia of hull structure about vertical axis through 
centroid of each section. 



K » ratio of average shear stress to maximum shear stress over 
a cross. -section for horizontal transverse shear. 



L 

M 

N 



- length overall of ship, ft. 



* bending moment, ft-lb. 
s an operational constant, 
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NOMENCLATUBE (Continued) 


p * 


an operational constant, * ^f~TT ~T7* 


Q « 


shear force, lb 


T » 


period of vibration, milliseconds 


w i ■ 


an operational constant * rn £ dx/L 


If 

£ 


an operational constant « J yf mo dx/L 


vv 


deflection patterns for a uniform prismatic bar obtained from 
consideration of simple bending theory only, according to Lord 
Kayleigh. 


f ■ 


frequency of vibration, cps 


g * 


. 2 

force of gravity, 32.2 ft/sec 


m * 


virtual mass per unit length, including structure, ballast, bee - 
flooding water, etc, , plus water entrained in horizontal vibration 
(ib sec 2 /* 2 ; 


« * 

X * 


mg , virtual weight per unit length (lb/ft) 

coordinate parallel to longitudinal axis of ship, measured from 
bow (ft) 


X * 


coordinate of longitudinal centroid of virtual mass 


y * 


coordinate parallel to beam of ship, measured from vertical 
plane through centroid of ship; measure of horizontal deflection 


y * 


the deflection pattern obtained from a cycle, before orthogonal 
lizatioa 


yO s 


value of unnormalized deflection curve taken at the normalizing 
point 
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NOMENCLATURE Continued 

? * coordinate parallel to vertical axis of ship measured from 

plane through centroid of ship. 

A a virtual displacement (pounds); see 

(9 * component of slope due to shear deflection 

T * component of slope due to bending 

? * component of slope due to bending before axis de terminating r+C 3 »TT 

yj ■ deflection before determination of axis*. * y 

r l' * slope before determination of axis; fj 4- * y* 

f « dummy variable, replacing x 

* circular frequency, (radians/ second) 

/ 

Note: Subscript to deflection indicates mode number . 

Hyphenated subscript indicates mode number and iterative cycle num- 
ber. 

Prime indicates differentiation with respect to x. 
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1. INTRODUCTION 

A. Historical 

For the last seventy years, the subject of vibrations of ship hulls 
has attracted ever increasing attention from the naval architect. Early 
investigators, notably Otto Schlick. emphasised the empirical treatment; 
but so complex is thedructure of a large ship, and so varied its design 
features, that it became clear that this approach would not iully serve the 
need. Ultimately then, we must look to the theoretical approach for the 
development of adequate formulations. 

Basically, the ship Is considered a slender free-free elastic beam, 
vibrating transversely to its long axis, for which the theory stems from 
Lord h&yleigh. f6l*From energy considerations he developed the basic 
equation based on simple bending. He then showed that for various 
assumed mode shapes, the one for which the energy balance gave the lowest 
frequency was the correct shape. He alao published a modified equation, 
accounting for the effect of rotatory inertia, but, considering this eLect 
negligible, did not emphasise it. 

Timoshenko formulated the shear flexibility term for JKayleigh's 
basic equation, and thin set down the well known Timoshenko mechanism, 
which considers the contribution of bending, shear, and rotatory inertia. 
While this formulation ia the most complete available, it is complex and 
difficult to solve exactly, even for a uniform beam, and in recent years a 
number of papers have neen devoted to this solution, generally involving 
Laplace transforms or scries solutions. 11, 33, 34, 63, et.al. 

Interest in the non-uniform beam inspired work along difterent lines, 

since a beam with a general distribution of mass, inertia, and section area 

is not amenable to exact solution. One approach is to attempt to for . ulatc 

one or more of these distributions into expressions suitable for analytic 
♦ Numbers in brackets refer to references listed in Appendix F. 
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treatment by operational calculus, e.g. £9, 10* 35* 37j . However* it 
ie evident that a reasonably accurate method of solution which is unre- 
stricted in its application, i. e. , can treat completely general distributions 
of structural characteristics, hss advantages over any isolated analytical 
solution* however rigorous. Thus we shall consider the iterative methods* 
both graphical and tabular. 

Stodola's work on turbines [8} suggested treating non-uniiorm beams 
by integrating Rayleigh's differential equation numerically or graphically. 
His methods also included the use of the funicular polygon* and Rayleigh's 
rule of minimum frequency. In using the successive integration method, 
he observed that while the two-noded mode behaved well, a three -noded 
starting assumption diverged during iteration until it returned to the two- 
noded mode; he failed to recognize the mutual orthogonality of normal 
modes. Viaoeilo Introduced this method into general use in the field of 
engineering, and J. L. Taylor and Schadlefsky applied It to ships (illustrated 
In £6fc] , as is Lewis's tabular method). All worked primarily with the 
first* or two-noded mode, and employed the simple -bending form of the 
equation of Rayleigh, though Taylor recognized the importance of Timo- 
shenko's shear term in ship vibrations. Koch C 2] proved the convergence 
of the method, and more important, demonstrated the necessity of 
rigorously maintaining orthogonality with all previous modes throughout 
the iterative method. 

Murray demonstrated the importance of transverse shear in ship 
vibrations £70], and introduced this effect into bcha dial sky's method, 
(rather than applying it as a subsequent correction to frequency as did 
McQpldrick £6&] ). _This refined method, first employed by Murray and 
the authors in work on a submerged submarine [71] , was noted also for 
a method of determining the slope axis by integration by parts of the 
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expression for orthogonality with rigid body rotation (pitch and yaw), and 
for calculation of higher modes than the first, facilitated by application of 
Koch's intermode orthogonality conditions. 

Other graphical methods have been employed by Southwell, Gumbel, 
Morrow, and others, but were suitable generally for the first mode only. 

Recently the principles employed in these iterative methods have ueen 
applied to lump parameter formulations for solution by sequence controlled 
digital calculators or by electrical analogy. £56, 67, 69] . la particular. 

Pro hi and Myklestad developed finite difference equations for use in the 
laylor Model Basin digital computer, taking account of simple bending, 
shear deflection, and rotatory inertia. Calculations made b ; this method 
for vibrution in the vertical plane for the attack transport U.S.5. NIAGARA 
with and without the secondary effects (shear flexibility and rotatory inertia) 
were compared with experimental data. I he increasing importance of 
secondary effects with increasing frequency was demonstrated, but for this 
vessel rotatory inertia was relatively unimportant In the vertical plane. 

B. Objective 

This investigation is concerned specifically with transverse vibration 
in the horizontal plane of a submerged submarine. Recent interest in the 
bodily response of a submarine to dynamic loading such as that accompanying 
a non-contact underwater explosion has led to the development of methods 
of modal analysis. These methods, and other demands of modern shi? 
design, have underscored the inadequacy of older established methods of 
determining natural modes of response of a ship. 

In particular, then, the object of this investigation was to establish 
and demonstrate a method of determining, as accurately as possible, the 
curvature pattern, deflection pattern, and natural frequencies of the normal 
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•nodes of interest in a transversely vibrating submarine of general mass, 
inertia, and section area distribution. Further, the method was to be 
capable of execution without resort to such expensive mechanisms as 
digital computers. 

C. Method: 

The method arrived at is basically similar to that of Schadlofsky 
f68^]> involving iterative cycles operating on an assumed deflection 
pattern, each cycle comprising the four successive intentions suggested 
by Rayleigh’s equation. The Timoshenko shear deflection correction is 
made to the slope o >tsned from bending, and the rotatory inertia correc- 
tion of Rayleigh is made to the derivative of the moment. Koch's ortho- 
gonality checks are made to filter out components of all lower modes, and 
Murray's slope axis determination is employed. Integration is performed 
by integraph (Plate 1), and calculations by commercial desk calculator. 

The method is demonstrated for the first three normal modes, using 
structural data from an existing three-eighths scale model of a submarine, 
including an allowance for entrained water. 

Shear lag. and the possible interaction between shear deflection 
and rotatory inertia are not considered. 

It is noted t at the problem is formulated as an eigenvalue prob- 
lem, in which solutions exist only for a discrete set of squared frequencies, 
each of which corresponds to a particular mode pattern, except that the 
zero eigenvalue indicates both translational and rotational vibration of 
the rigid body. 
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II PROCEDUEi, 



The basic theory of vibration of a non-uniform, transversely 
vibrating beam as formulated by L-ord R? ylelgh, considering only the 
simple bending effect, can be expressed briefly by: 

'b* ( £1 tly ) « - (1) 

where the assumption of simple harmonic motion ot circular frequency 
has been incorporated. Four successive integrations with respect to 
this equation comprise the framework upon which the iterative method of 
this investigation has been developed, ’x he deflection and frequency 
included in the inertia loading term on the right side of (1) are assumed 
values, and the deflection indicated on the left, obtained by the quadruple 
integration, is the first estimate of the true deflection pattern. Ihe first 
estimate of the square of the true circular frequency is obtained uy 
comparison of the assumed and derived deflection antipodal amplitudes, i.e.. 



( *>* ) 



final 



* S 63*y*) star t 
(y*) final 



( 2 ) 



Note also that the second integral of inertia loading, viicb is the curve of 
bending moment, must be divided by the curve ©. local flexural rigidity 
before the last two integrations are performed, since Li is a variable in x. 

We first modify the process by applyiug to the shear curve yielded by 
the first integration a correction for rotatory inertia, as follows: 



dM a Q - _Imr (3) 

dx A 

where Y is the bending component of slope, bince this quantity occurs 
later in any iterative cycle than M, we must use the T value from the 
previous cycle as a best approximation, ’thus this correction is not made 
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until the second cycle of each mode. 

The second modification of the process is the correction of the slope 
curve obtained by integrating the curvature distribution. This correction 
introduces the effect of shear flexibility, and the total slope resulting is 
expressed by: 



Here the shear vsiia Q has been obtained before the slope calculation, 
and the correction can be made in each cycle. 

Application of boundary conditions and orthogonality conditions, 
including orthogonality (with respect to mass) with the aero eigenvalue 
motions and with all modes below the one being obtained, at the 
appropriate points in and after each cycle, determine the four constants 
of integration, and refine the pattern obtained iron-: the iterative cycle. 

In this investigation, since mass and flexural and shear rigidities 
have general, hence irregular distributions, analytical integration was 
not possible. The process was therefore carried out in each instance by 
mans of * mechanical integraph with a 23" carriage beam (manufactured 
in Switzerland by C. Coradi). 1 his instrument draws sn integral curve 
when its stylus is made to trace the argument curve (Plate I). All calcu- 
lations have been carried out and plotted at intervals of one-fortieth of 
ship length, as this was found to be the greatest interval which would 
detect important variations in the functions used. 

The sequence of operations is aafollows: 

(1) 1 he structural characteristics needed are computed fron. 

Available plana and tabulated (Appendix E). 




( 4 ) 



where 



e - q 



( 3 ) 
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(2) A tabulation of functions of structural characteristics and of 
the centroid of virtual mass which will be used throughout 
the work is prepared. (Table IV, Appendix C). To this 
table are added a tabulation of a function of each modal 
deflection pattern obtained, for use in orthogoaulizing sub- 
sequent modes. 

(3) A tabulation of constants computed for repeated use is prepared 
(Table V, Appendix C). 

(4) A tabulation of the Bayieigh mode patterns (obtained from 
application of his simple -bending theory for a uniform slender 
beam) is prepared, for all modes to be obtained. These are 
normalized patterns. 

(5) The Bayieigh pattern for the mode ^eiag calculated is orthogona- 
Lized with translation, rotation, and all previous modes (there 
would be none for the two-noded mode). In practice, the Bayieigh 
pattern need not be used, any curve of the correct general snape 
is suitable for a starting guess, and the better the guess, the more 
rapid the conversion to an acceptable solution. Perhaps the best 
choice would be the pattern freer, a similar ship, it such is 
available. 

The resulting orthogonalized pattern will now be used as a 
starting point of the iterative process. 

(6) Inertia loading based on this pattern is computed and plottel** a 
function of X/L along the length of the ship. 

(7) This is integrated to obtain shear, to which is added the rotatory 
inertia term, giving the derivative of bending moment. (In the 
first cycle of each mode shear is the differential of bending 
moment). M' is plotted. 
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(3) M 1 is integrated, and the integral curve is divided, at each 
station, by flexural rigidity to obtain curvature, which ia 
plotted. 

(9) Curvature is integrated to obtain the bending component of slope 
to within a constant. 

(10) Shear is divided by shear rigidity, and the resulting shear slope 
is added to the bending cc uipaucat, and the total is used in an 
auxiliary process to determine a constant which locates the 
axis (the constant of the integration in (9)), which is added also, 
to give total slope . This is plotted. 

(11) Slo^e is integrated to give the deflection pattern to within an 
additive constant, which is used in an auxiliary process to obtain 
the axis shift. The resultant deflection pattern is tabulated. 

(12) This pattern is now submitted to the same orthogonality checks 
used on the starting guess. I he corrected curve is normalised 
to give a pattern whose value is unity at the center - or center' 
most -alter -antinode. The unnormalized value at this antioode 
is retained for use in the frequency calculation. The normalised 
pattern becomes the oasis of the inertia loading for the next 
cycle, or is accented as the final pattern. This decision is based 
on the amount of change of deflection pattern and frequency be- 
tween this cycle and the previous one. 

(13) Frequency is obtained by equation (2), where ( y*)#tart a 1» aud 
(y*) final is the unaormalized antinode deflection retained in (12). 

(14) In the final cycle of each mode, the curvature pattern is 
normalized by dividing throughout by (y*) final and retained as 
the final curvature pattern. 

(15) The shear slope distribution is subtracted from the corrected 
total slope distribution to give the corrected bending slope. This 
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is normalized bydivision by (y*) final and retained for 
use in the rotatory inertia correction of the subsequent 
cycle. 

(16) The pattern obtained in (12) is now used to repeat the 
entire cycle, if deemed necessary. 

Appendix A, Part 1 contains the derivations of the important relations 
used. Part 2 is a detailed outline of the procedure, containing formulations, 
of all operations as used in the application of the method illustrated in 
Appendices B and C. 

Ihe method was executed using structural data for a three -eighths 
scale model oi a submarine, and three lateral modes were computed, each 
having been carried through three cycles of iterations, this being sufficient 
to show convergence of the method and to reduce residual error, as 
indicated roughlyby the last increment of change, to an acceptable level. 
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111 Kt.5 VLT 

by observation of ike convergence of the results of iteration, 
patterns, normalized to give a value of unity at the center or 
centt r most after anti node, are presented in Figure XXXIIi. The 
v or responding curvature patterns, front which strain patterns can 
be obtained, were normalised in such a manner that when the factor 
between normalised deflection and true deflection maxima has been 
obtained, this same factor multiplied by the normalised curvature pattern 
will give the maxima of the true curvature pattern, all dimensions 
being consistent. These patterns are presented in Figure XXXIV. 

Table I gives the final frequencies and periods of the three modes 
calculated for the model: 

TAj&Lk I 

MATUEAu TKANSVfcfcSC FFJtQUFNCILS AND PSEIQDS 
FOh IUO. 5 FOOT SUBMXKGhD SUbM^&INL MODbL 



Mode 


Natural F xeq. (cpm) 


Katio 


Period (milliseconds) 


1 


417 


1 


143.9 


2 


S3* 


2.01 


7i. 5 


3 


1391 * 


3. 34 


43. 1 



* See supplementary discussion. Appendix D 

Table II lists the offsets of the normalized deflection and curvature 
patterns pictured in Figures XXXIII and XXXIV. 

to the method developed, it was found to be convergent, though 
it appears that the delayed introduction of rotatory inertia into the pro 
ce»t must act to delay convergence. This is compensated for by the 
frequent orthogonality checks, especially their application to th^ 
starting pattern. 

The results and method are discussed in the subsequent section. 



5 * 



J Ul 



,«el> . . : I 1 • “ • "•!'* • - • • i ‘* i •'***'* X 

.. . • , j ]4 • .. 4 ' * ' 1 • h ‘ * * ,ln 0 ' 1 

• «vJ . 1 ./ - »♦*»*'- * *H. ... 

, *)' »• c <-f «* 1 *■ • •• ■••••» *•■* • 

» 4 ,.i •* *• , J..I* ' • ’ “• W. .x • * •• • * i *' 1 *- 1 * 

Jk , , sm V *11 iit • * **'♦ l •** • s 

. '••»- -• ' •** *• ' 

,.*» J . u i • ♦•»' * - • »* -A'** • l» • ,4 ‘ ' ■ ^ 

riA.A • . . •* fcM* •**.•! * •** 

k 4 rt,. j « .♦! t « V »; W«» • • * - * * * ‘ * 1 

i »>*#<%■ •!» i . 



C .i% 4 ‘ .K^’J * - • ' M • 1 

. *c»<. ... . .-l . ► . *v c -t i» "-*!••• 1 • • * 



lit ) fcillK 


• it» t 


I*** 


•wu 


f.Ul 




Til 


i 


• JV 


A.l 


hi 


V 


1 .!• 


M .t 


* »♦<-• 


1 




a i 


^ ,• • ,J 


o .. V »• • 



' ' ‘ .Oft ••' l< ' “* ^ **" 

, . . y 11. ..*•'•• • i • «* .41 4. ** •**• ‘ 



• « * 



, < -a-r Vi «#>• -'*•'* • : -** /-.*-» ••» 

., t '> * • i<» 1 . * "* • 1 » b * l * *** 

,* ,.r.H ,*• -A. Vi * '«*•** 

.r • : »». i *•» »* «=>• 



• 3 « f •• n * . . 



MV 



,•. ,i »o • • *•• • • • i 



1 ■» ’ * 



22 



OHDIHATI^ or pr.FLSCTXQM AfO CUaVA?UriS PATTEHna 



x/t 


*1 


1 i 

* 1 




n 

' 2 


f 3 


1 1 

* 3 


o.ooo 


-3.408 


0 


+3.041 


0 


+3*925 


0 


0*025 


3.054 


- 258 


2.578 


+ 349 


3*116 


+ 1916 


0.050 


2.695 


w 


2.101 


2440 


2.231 


9581 


0.075 


2.347 


1614 


1.609 


4792 


1.630 


14371 


0.100 


1.994 


1840 


1.161 


6142 


0.729 


15696 


0.125 


1.641 


2117 


0.726 


6515 


+0.102 


16968 


0.150 


1.299 


2363 


+0.336 


6 ?6l 


-0.417 


16727 


0.175 


0.988 


2554 


-0.018 


6937 


0.808 


15751 


0.200 


0.682 


2311 


0.307 


7316 


1.074 


15005 


0.225 


0.395 


3076 


0.567 


7506 


1.178 


13877 


0.250 


-0.125 


3374 


0.762 


7589 


1.150 


11536 


0.275 


+0.10> 


3644 


0.891 


7527 


1.071 


9337 


0.300 


0.325 


4054 


O.963 


7541 


O.662 


6815 


0.325 

0.350 


0.511 

0.561 


im 


0.977 

0.940 


7474 

7074 


0.596 
-0, 295 


3S96 
+ 489 


0.375 


0,788 




0.845 


6138 


*0.036 


- 2801 


0.400 


0.378 


6*468 


0.706 


4352 


0.330 


*4552 


0.425 


0.944 


3240 


0.532 


2295 


0.602 


4899 


0.450 


0.987 


2336 


0.336 


1317 


0*793 


3445 


0.475 


1.000 


2032 


-0.133 


* 588 


0*928 


5977 


O.5OO 


0.989 


2042 


+0.005 


- 154 


0.996 


6X30 


0.525 


0.965 


3361 


O.302 


1076 


0.975 


6860 


0.550 


0.912 


451? 


0.519 


2627 


0.875 


8160 


0.575 


O.835 


5X31 


0.693 


4400 


0.677 


7271 


0.600 


0.711 


5019 


0.833 


5625 


0.617 


4458 


0.625 


O.568 


4763 


0.938 


65 76 


+0,083 


— 86l 


0.650 


0.383 


4502 


0.995 


7397 


-0.275 


+ 3024 


0.675 


+0.169 


430r 


0.98? 


8124 


0.655 


6922 


0*700 


-9^087 


4010 


0.920 


8653 


1.002 


10*36 ' 


0.725 


0.350 


3727 


0.780 


8360 


1.227 


14653 


0.750 


O.654 


3367 


O.502 


8808 


1*373 


17691 


0.775 


0.976 


3045 


+Q.304 


8570 


1*389 


20157 


0.800 


1.310 




—0#02*s 


7845 


1.281 


21278 


0.025 


1.662 


2365 


0.411 


7405 


0.951 


22813 


0.050 


2.032 


2259 


0.849 


7406 


-0,687 


25419 


0.875 

0.900 


2.389 

2.783 


2773 

3341 


1.337 

1.883 


9187 

10760 


+0.155 

1.044 


34836 

44791 


0.925 


3.207 


2360 


2*503 


7843 


2.297 


39530 


0.950 


3.624 


1036 


3*174 


2091 


3*711 


21556 


0.975 


4.066 


- 129 


3.803 


-5229 


5*154 


+ 1437 


1.000 


-4.502 


0 


-4.450 


0 


+6. 590 


0 









' ,f - 
' » 

- # 

*• * , r 

• i.r 

. 

>r<r 4 

i t i 4 '-*4 




>v * 




401 



25 



IV. DISCUSSION OF RESULTS 
A . Mode Patterns: 

The shape of the mode pattern* (Figure XXXIIi) generally n tt the 
authof*' expectations, in ships and submarines, uoth the mass and 
fiexu al rigidity are greater at the midsection than at the ends. This 
has the effect of altering the mode shape so that the center deflections 
are smaller relative to the end deflections than in a uniform beam, 
and the nodal points are closer to the midship section. For example, 
for the first mode, the center deflection maximum of the uniform beam 
amounts to 60% of the end deflection, while in the case of the submarine 
studied, the anti node deflection was 29.47* of the bow deflection and 
22.8% of the stern deflection. The nodes for the uniform beam ocvur 
at 22% i rorn the ends, while for the submarine investigated, they 
occurred at 26.2% from the bow and 30. 7% from the stem respectively. 

Figure I clearly demonstrates the pronajxed increase in mass, 
sectional area, moment of inertia, and shear factor K in the vlvinity of 
the hard tank and conning tower amidships in the case of the submarine. 
One can readily arrive at two consequences of this discontinuity of 
structural strength. The first of these has already been mentioned - 
that of the very high end deflections relative to those in the middle body. 
This behavior will obviously occur for ail the modes, with the resulting 
trend towards patterns with many small loops in the middle with large 
end branches. It is easy to see how any loading, such as an explosion, 
occurring in the middle body of a submarine can incur very large deflec- 
tions at the ends with attendant small deflection near the load. In the 
case of the third mode, the deflection at the stern will be six times that 

at the center I The reader will readily appreciate the importance 
this characteristic of the mode patterns with respect to the bodily 
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response of submarines. Severe whipping may be expected to occur 
whenever the contributions of the odd and/or even numbered mode end 
deflections become additive for any general vibratory motion of the 
submarine . 

The second consequence of the irregular distribution of structural 
characteristics becomes manifest in the effect of rotatory inettia and 
shear flexibility. For a given structure mass per unit length, a sub- 
marine will be many times stiffer than a solid prismatic bar, since 
in the case of the former, the crossection is essentially a circular ring 
of material relatively distant from its centroid. This fact, coupled with 
that of the nature of the longitudinal distribution of flexural rigidity, 
explains why the magnitude of the rotatory inertia correction to the 
frequency of the first mode is somewhat greater for the submarine than 
for the uniform beam or even a surface ship. Its increased moment of 
inertia for a given mass, ft is compared to the uniform beam results in 

larger values of Im . Even for equal slopes, the rotatory inertia 

A 

correction would be larger for the submarine. However, the in- 
creased flexibility of the end sections over their uniform beam counter- 
parts causes the slope to increase over all the length of the submarine, 
the greatest increase occurring between the ends and the nodes. The 
steepness of all the deflection patterns near the ends results in a rather 
selective effect of rotatory inertia upon the patterns. In general, the 
effect of rotatory inertia upon the mode patterns studied was most pro- 
nounced at the ends, and resulted in an increase in the relative deflection* 
A comparison of experimental and theoretical mode patterns for a similar 
submarine where rotatory inertia was neglected in the latter calculation, 
shows larger experimental end deflections than the theory predicted. 

The authors feel that such a discrepancy may be due to the application 
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